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Abstract 

We provide a snapshot of Dyson-Schwinger equation applications to the the- 
ory and phenomenology of hadrons. Exact results for pseudoscalar mesons 
are highlighted, with details relating to the Ua{X) problem. Calculated 
masses of the lightest J = 0, 1 states are discussed. We recapitulate upon 
studies of nucleon properties and give a perspective on the contribution of 
quark orbital angular momentum to the spin of a nucleon at rest. 

1 Introduction 

Numerous salient features exhibited by the physics of mesons and nucleons 
arise nonperturbatively in QCD. Two phenomena strike one immediately: 
confinement and dynamical chiral symmetry breaking (DCSB). DCSB is the 
better understood of these emergent phenomena; e.g., it explains the origin 
of constituent-quark masses and underlies the success of chiral effective field 
theory. Confinement, on the other hand, remains only an empirical fact; 
viz., colored objects have not hitherto been observed in isolation. A fact too 
often ignored is that the potential between infinitely heavy quarks measured 
in numerical simulations of quenched lattice-QCD - the static potential - is 
not related in any known way to light-quark confinement. 



2 DCSB 

Understanding DCSB within QCD proceeds from the gap equation [1]; namely, 
the Dyson-Schwinger equation (DSE) for the dressed-quark propagator: 

S{p)-' = Z,it^-p + m^-) + Z, pg^D^^ip - q)^^^S{q)T:{q,p), (1) 
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where represents a Poincare invariant regularization of the integral, with 
A the regularization mass-scale, D^,, is the dressed-gluon propagator, Vy is 
the dressed-quark-gluon vertex, and m*^™ is the quark's A-dependent bare 
current-mass. The vertex and quark wave-function renormalization con- 
stants, Zi 2(C^5A^), depend on the gauge parameter. 
The solution of Eq. ([T]) can be written: 

Siv) = ' ^^^''^'^ (2) 

wherein ('^) is the wave-function renormalization and M(p^) is the mass 
function. The latter is independent of the renormalization point, (. The 
dressed propagator is obtained from Eq. ([T]) augmented by the renormaliza- 
tion condition 5" 1^2=^2 = ij ■ p + "^(C^) , where m{(^) is the running 
mass: Z2(C^, A^) m^™(A) = Z4(C^, A^) m(C^) , with Z4 the Lagrangian-mass 
renormalization constant. In QCD the chiral limit is strictly defined by [1]: 
Z2(C^, A^) m^™(A) = 0,VA^ ^ C^, which states that the renormalization- 
point-invariant current-quark mass m = 0. 

Only in the chiral limit is it possible to unambiguously define the gauge 
invariant vacuum quark condensate in terms of S{p) [1-3]. This emphasizes 
that gauge covariant quantities contain gauge invariant information. The 
condensate is the order parameter most commonly cited in connection with 
DCSB. Nonetheless, M(p^) is a more fundamental indicator: the condensate 
is only a small part of the information it contains. 

In perturbation theory it is impossible in the chiral limit to obtain M(p^) 7^ 
0: the generation of mass from nothing is an essentially nonperturbative 
phenomenon. On the other hand, it is a longstanding prediction of non- 
perturbative DSE studies that DCSB will occur so long as the integrated 
infrared strength possessed by the gap equation's kernel exceeds some criti- 
cal value [4]. There are strong indications that this condition is satisfied in 
QCD [5,6]. It follows that the quark-parton of QCD acquires a momentum- 
dependent mass function, which at infrared momenta is ~ 100-times larger 
than the current-quark mass. This effect owes primarily to a dense cloud of 
gluons that clothes a low- momentum quark [7]. 

A great deal has been learnt from the gap equation alone. To highlight 
only one recent example [3], realistic Ansdtze for the the gap equation's kernel 
indicate that there is a critical current-quark mass, rhcr, above which M{p'^) 
does not possess an expansion in m around its chiral-limit value. For a pion- 
like meson constituted from a quark, /, with mass rhcr and an equal- mass 
different-flavor antiquark, ^, m^^ = 0.45 GeV. Since physical observables, 
such as the leptonic decay constant, are expressed in terms of M(p^), it 
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follows that a chiral expansion is meaningful only for (m^^)^ < 0.2 GeV^. 
This entails, e.g., that it is only valid to employ chiral perturbation theory to 
fit and extrapolate results from lattice-regularized QCD when the simulation 
parameters provide for < 0.2 GeV^. Lattice results at larger pion masses 
are not within the domain of convergence of chiral perturbation theory. 

3 Mesons 

For 0~ mesons the axial-vector Ward-Takahashi identity is of fundamental 
importance, and recently its implications for neutral pseudoscalars and rj-rj' 
mixing have been elucidated [8]. In the general case the identity is written 

P,n^{k- P) = 5-i(fc+)^75^'^ + il.T'^S-^k.) - 2iM'''Tl{k- P) - A\k- P) , 

(3) 

wherein: {J^°'\ a = 0, . . . , N'j — 1} are the generators of U{Nf); the dressed- 
quark propagator S = diag[S'u, Sd, Ss, Sc, Sb, . . .]; A1(C) is the matrix of run- 
ning current-quark masses and Ad""^ = tip [{J-'"', , where the trace 
is over flavour indices. The inhomogeneous axial- vector vertex in Eq. ([3]), 
^tfj-i^'i -^)' where P is the total momentum of the quark-antiquark pair and k 
the relative momentum, satisfies a Bethe-Salpeter equation (BSE), and hke- 
wise the pseudoscalar vertex, T\{k;P). The final term in Eq. ([3]) expresses 
the axial anomaly. It involves 

Au{k;P)=Jd^xd'ye'^'+-^-''-y'>Nf{j^^q{x)Q{0)q{y)), (4) 

wherein the matrix element represents an operator expectation value in full 
QCD and Q{x) = €.^,yp„F^^F^^{x) = d^K^{x) is the topological charge 
density operator, where F^^, is the gluon field strength tensor|l| 

In considering the ?7A(l)-problem one need only focus on the case A'^ ^ 
because if that is false, then following Ref. [1] it is clear that the t]' is certainly 
a Goldstone mode. A^ is a pseudoscalar vertex and can therefore be expressed 

A\k-P) = j^'j,[zSAk;P) + i-P:FA{k;P) 

+^-kk- Pgj,{k- P) + a^^k^PMAik; P)] . (5) 

Equation can now be used to derive a collection of chiral-limit, point- 
wise Goldberger-Treiman relations, important amongst which is the identity 

^NB. While Q{x) is gauge invariant, the associated Chern-Simons current, iiT^, is not. 
Thus in QCD no physical state can couple to iC^. Hence, physical states cannot provide a 
resolution of the so-called JTa (l)-problem; namely, they cannot play any role in ensuring 
that the if is not a Goldstone mode. 
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Table 1: Masses (GeV) of the lightest J = 0, 1 states produced by the 
rainbow-ladder DSE truncation of Refs. [17, 18] with the parameter values: 
u) = 0.4 GeV, ujD = (0.72 GeV)^; and current-quark masses m„^rf(l GeV) = 
5.45 MeV, ms(lGeV) = 125 MeV. The rainbow-ladder kernel gives ideal 
flavour mixing for all states. See the text for further discussion. 
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2fn'En'{k; 0) = 2i?o(/c^) — £A{k] 0) , where Bo^k"^) is obtained in solving the 
chiral-limit gap equation. It is plain that if 

£:^(fc;0) = 25o(A;2), (6) 

then fri'Eri'{k; 0) = 0. This being true, then the homogeneous Bethe-Salpeter 
equation for the t]' does not possess a massless solution in the chiral limit. 
The converse is also true. Hence Eq. ([H]) is a necessary and sufficient condi- 
tion for the absence of a massless r]' bound-state. The chiral limit is being 
discussed, in which case Bo{k'^) ^ if, and only if, chiral symmetry is dynam- 
ically broken. Thus the absence of a massless rj' bound-state is only assured 
through the existence of an intimate connection between DCSB and an ex- 
pectation value of the topological charge density. A relationship between the 
mechanism underlying DCSB and the absence of a ninth Goldstone boson 
was also discussed in Ref. [9]. 

Reference [8] also derives corollaries, amongst which are mass formulae 
for neutral pseudoscalars, and presents an Ansatz for the Bethe-Salpeter 
kernel that enables their illustration. The model is elucidative and phe- 
nomenologically efficacious; e.g., it predicts tj-t]' mixing angles of ~ —15° 
and 71^ -r] angles of ~ 1°; and suggests a strong neutron-proton mass differ- 
ence of 0.75 {nid — rriu). 

The use of DSEs to study meson phenomena is empowered by the exis- 
tence of a systematic, nonperturbative and symmetry-preserving truncation 
scheme [10,11]. It means that exact results, such as those indicated above, 
and others related to radial excitations and/or hybrids [12-14], and heavy- 
light [15] and heavy-heavy mesons [16], can be proved and illustrated. 

In the latter connection, the renormalization-group-improved rainbow- 
ladder truncation of the gap and Bethe-Salpeter equations introduced in 
Refs. [17, 18] has been widely employed. To exemplify that, in Table[T] we 
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report calculated results for the masses of the lightest J = 0, 1 states [19]. It 
is true in general that the truncation is accurate for the and 1 light- 
quark meson ground states. In these channels it can be seen algebraically 
that contributions beyond rainbow-ladder largely cancel between themselves 
owing to Eq. ([S]) [11,20-22]. The remaining columns in the table deserve 
special attention because they show clearly the path toward improvement. 

Terms beyond the rainbow-ladder truncation are known to add construc- 
tively in the 0"'""'" channel [23]. Hence the leading order truncation is a priori 
not expected to provide a good approximation. Further understanding is 
provided by an exploration of the contribution from two-pion intermediate 
states to the mass and width of this lowest-mass scalar. A rudimentary anal- 
ysis shows that a realistic description is attainable therewith [25]; viz., it 
gives a pole position = (0.578 — i 0.311) GeV. This is not the end of 
the scalar story but it is a sensible path to follow, in particular because a 
QCD-level mechanism is precisely specifiedH 

Compared with experiment, the masses of the axial- vector mesons l"*"^ are 
poorly described by the rainbow-ladder truncation: ~ 400 MeV of repulsion 
is missing from the kernel. A cruder model does better [26,27]. The latter 
studies and a more recent analysis [28] indicate that at least part of the 
defect owes to the absence of spin-flip contributions at leading-order. Such 
contributions appear at all higher orders and are enhanced by the strongly 
dressed quark mass function. It is in this way that the meson spectrum 
can be used to probe the long-range part of the light-quark interaction and 
thereby to chart the nonperturbative behavior of QCD's /9-function. 

The last three columns describe systems with so-called exotic quantum 
numbers. Of course, these states are exotic only in the context of the naive 
constituent quark model. In QCD they correspond simply to interpolating 
fields with some gluon content and are easily accessible via the BSE [29]. 
Nonetheless, while the rainbow-ladder truncation binds in these channels, 
the shortcomings encountered in the l"*" channels are also evident here, for 
much the same reasons. Reliable predictions for the masses of such states 
will only be obtained once improved kernels are developed. At the very least, 
one must have reliable predictions for axial-vector masses before drawing any 
conclusions about the so-called exotics. 

One might pose the question of whether, in the context of bound-state 
studies in which model assumptions are made regarding the nature of the 
long-range interaction between light quarks, anything is gained by working 
solely with Schwinger functions. This means, in part, constraining oneself to 

^In rainbow-ladder truncation, at least up to the c-quark mass, the ordering of meson 
masses is 0^+(lS') < 0++(lS') < 0-+{2S) < 0++{2S) [24]. 
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work only with information obtained from the DSEs at spacelike momentajfl 
According to a recent study [30] the answer is no. It analysed the capacity 
of Schwinger functions to yield information about bound states, and estab- 
lished that for the ground state in a given channel the mass and residue are 
accessible via rudimentary methods. However, simple methods cannot pro- 
vide dependable information about more massive states in a given channel. 
Indeed, there is no easy way to extract such information. An approach based 
on a correlator matrix can be successful but only if the operators are care- 
fully constructed so as to have large overlap with states of interest in a given 
channel and statistical errors can be made small; viz., ~ 1%. While it is pos- 
sible in principle to satisfy these constraints, doing so is labor intensive and 
time consuming. That is only justified in the absence of model-dependence. 

4 Nucleons 

The discussion of DSEs at MENU04 did not describe the study of nucleons, 
stating only that it was feasible [31]. Material progress has been made in the 
interim. We now possess a level of expertise roughly equivalent to that we had 
with mesons approximately ten years ago; viz., phenomenology constrained 
by the significant body of knowledge we have gained in meson applications. 

The nucleon appears as a pole in a six-point quark Green function. The 
pole's residue is proportional to the nucleon's Faddeev amplitude, which 
is obtained from a Poincare covariant Faddeev equation that adds-up all 
possible quantum field theoretical exchanges and interactions that can take 
place between three dressed-quarks. Poincare covariance is crucial because 
modern experimental facilities employ large momentum transfer reactions. 

A tractable truncation of the Faddeev equation is based [32] on the ob- 
servation that an interaction which describes mesons also generates colour-3 
diquark correlations [33]. For ground state octet and decuplet baryons the 
dominant correlations are 0"*" and l"*" diquarks because, e.g.: the associated 
mass-scales are smaller than the baryons' masses [26,34,35], namely (in GeV) 
"^Mo+ = ~ 0-^' = = ^{dd)^+ = 0.9 - 1.0; and the elec- 

tromagnetic size of these correlations is less than that of the proton [36] - 
^Mo+ ^ 0.7 fm, which implies ~ 0.8 fm based on the p-meson/vr- 

meson radius-ratio [37,38]. 

The Faddeev equation's kernel is completed by specifying that the quarks 

■^Lattice-regularized QCD provides a background to this question. That approach is 
grounded on the Euchdean space functional integral. Schwinger functions; i.e., propagators 
and vertices at spacelike momenta, are all that it can directly provide. It can only be useful 
if methods are found so that the question can be answered in the affirmative. 
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are dressed, with two of the three dressed-quarks correlated always as a 
colour-3 diquark. Binding is then effected by the iterated exchange of roles 
between the bystander and diquark-participant quarks. A Ward-Takahashi- 
identity-preserving electromagnetic current for the baryon thus constituted 
is subsequently derived [39]. It depends on the electromagnetic properties of 
the axial-vector diquark correlation. 

A study of the nucleon's mass and the effect on this of a pseudoscalar me- 
son cloud are detailed in [40] . Lessons learnt were employed in a series of stud- 
ies of nucleon properties, including form factors [41-44]. The calculated ratio 
/ipG'^(g2)/G'P^(Q2) pg^gggg through zero at ~ g.S GeV^ [42]. For the neu- 

2 

tron, in the neighbourhood of Q"^ = 0, /i„ G'^{Q'^)/Gl.j{Q'^) = Q^, where 
r„ is the neutron's electric radius [43]. The evolution of UpG^^iQ"^) / G^j^^^Q"^) 
and lJ^nG'^{Q'^) / G1^{Q'^) on Q"^ > 2 GeV^ are both primarily determined 
by the quark-core of the nucleon. While the proton ratio decreases uni- 
formly on this domain [41,42], the neutron ratio increases steadily until 
Q'^ ~ 8 GeV^ [43]. A comparison of the Pauli/Dirac form factor ratios for 
the neutron and proton is presented in Ref. [7]. 

Of significant interest is the distribution of an hadron's spin over the 
quark constituents and their angular momentum. In a Poincare covariant 
approach that can be calculated in any frame. The rest frame is physically 
most natural. The pion was considered in Ref. [16]. The answer is more com- 
plicated for the spin-i nucleon. In the truncation just described a nucleon's 
Faddeev wave-function is expressed through eight scalar functions: no more 
are needed, no number fewer is complete. Two are associated with the 0^ 
diquark correlation: iSi 2, and six with the l"*" correlation: Ai^___fi. In the rest 
frame in this basis one can derive the following "good" angular momentum 
and spin assignments, which add vectorially to give a J = ^ nucleonQ 



L = 0,S=l 


L = l,S=l 


L = 1,S=I 


L = 2,5 = 1 




S2 ,Ai ,82 , 


C2 





(7) 



Bl = + , B2 = + 

C, = C2 = A^-Aq. 



These assignments are straightforward to demonstrate and understand; e.g, 
in the rest frame of a relativistic constituent quark model the iSi,2 terms 
correspond, respectively, to the upper and lower components of the nucleon's 
spinor. 

To exhibit the importance of the various L-S correlations within the nu- 
cleon's Faddeev wave-function we report the breakdown of contributions to 



^Equation (3.35) or Ref. [45] contradicts Fig. 6 of that reference. Equation ^ herein 
describes the correct assignments. 
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the nucleoli's canonical normalization 
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To illustrate how to read Eq. Q we note that the largest single entry is as- 
sociated with Si® Si, which represents the quark outside the scalar diquark 
correlation carrying all the nucleon's spin. That is the w-quark in the proton. 
However, it is noteworthy that a contribution of similar magnitude is associ- 
ated with the axial- vector diquark correlations, expressing mixing between p- 
and (i- waves; viz., Ci ® C2 + C2 (S> Ci. With C2 all quark spins are aligned with 
that of the nucleon and the unit of angular momentum is opposed, while with 
Ci all quark spins are opposed and the two units of angular momentum are 
aligned. This contribution is more important than those associated with ^2; 
namely, scalar diquark terms with the bystander quark's spin antiparallel. 
Finally, for the present, in this context one single number is perhaps most 
telling: the contribution to the normalization from (L = 0) (i^ = 0) terms 
is only 37% of the total. 



5 Coda 

The DSEs provide a natural vehicle for the exploration of confinement and 
DCSB. DCSB is a remarkably effective mass generating mechanism. For 
light-quarks it is far more important than the Higgs mechanism. It is under- 
stood via QCD's gap equation, which delivers a quark mass function with a 
momentum-dependence that connects the perturbative domain with the non- 
perturbative, constituent- quark domain. The existence of a sensible trunca- 
tion scheme enables the proof of exact results using the DSEs. The scheme 
is also tractable. Hence the results can be illustrated and predictions made 
for observables. The consequent opportunities for rapid feedback between 

^The entry in location Si ® Si indicates the integrated contribution associated with 
iSf . The entries are reweighted such that the sum of the squares of the entries equals one. 
Positions without an entry are zero to two decimal places. 
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experiment and theory brings within reach an intuitive understanding of 
nonperturbative strong interaction phenomena. 

It can be argued that confinement is expressed in the analyticity proper- 
ties of dressed Schwinger functions [4] . To build understanding it is essential 
to work toward an accurate map of the confinement force between light- 
quarks. Among the rewards are a clear connection between confinement 
and DCSB, an accounting of the distribution of mass within hadrons, and a 
realistic picture of hybrids and exotics. 

It is important to understand the relationship between parton properties 
on the light-front and the rest frame structure of hadrons. This is a challenge 
because, e.g., DCSB, a keystone of low-energy QCD, has not been realized in 
the light-front formulation. Parton distribution functions must be calculated 
in order to learn their content. Parametrization is insufficient. It would 
be very interesting to know how, if at all, the distribution functions of a 
Goldstone mode differ from those of other hadrons. Answers to these and 
kindred questions are being sought using the DSEs [46,47]. 
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